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I. INTRODUCTION 

The asymptotic elimination of fast time scales from the 
Vlasov equation presents important analytical and com- 
putational advantages for its solution in complex plasma 
geometries 0. Here, fast time scales in a collisionless 
plasma are either associated vifith particle orbital dy- 
namics (e.g., the fast gyromotion of a charged particle 
about a magnetic field line) or wave-particle dynamics 
(e.g., the fast oscillatory motion of a charged particle in 
the presence of a high- frequency electromagnetic wave). 
In the former case, the elimination of a fast orbital time 
scale is often associated with the construction of an adi- 
abatic (action) invariant (e.g., the magnetic moment of 
a charged particle in a strong magnetic field). Next, the 
construction of a reduced Vlasov equation, from which 
one or more fast time scales have been removed, can 
either be obtained through an iterative solution of a 
perturbatively-expanded Vlasov equation j2] or by per- 
forming one or more near-identity phase-space transfor- 
mations resulting from applications of Hamiltonian per- 
turbation theory 3] . The present paper focuses on appli- 
cations of Lie-transform Hamiltonian perturbation the- 
ory. 

The most general setting for carrying out Hamilto- 
nian perturbation theory 4] on the Vlasov-Maxwell equa- 
tions is to use an eight-dimensional extended phase space 
with coordinates Z°- ~ (z; w, t), where z denotes regular 
(six-dimensional) phase-space coordinates, w denotes the 
energy coordinate, and time t denotes its canonically- 
conjugate coordinate. The extended Hamilton's equa- 
tions (r denotes the Hamiltonian-orbit parameter) 



= {z^ n} = r\z] 



dnjz) 

dZ^ 



(1) 



are expressed in extended phase space in terms of an 
extended Hamiltonian Ti{Z) = H{z, t) — w, where 
H{z,t) denotes the regular Hamiltonian, and the ex- 
tended phase-space Lagrangian F = Fa dZ", which is 
a differential one-form in extended phase space (summa- 
tion over repeated indices is, henceforth, implied). Note 



that the physical motion in extended phase space takes 
place on the subspace H = 0: 



H{z,t). 



(2) 



The extended Poisson bracket { , } is obtained from the 
extended phase-space Lagrangian F, first, by construct- 
ing the Lagrange matrix w (with antisymmetric compo- 
nents Wab = da^i, — dbTa) associatcd with the differen- 
tial two-form ti) = dF = i Uab dZ" A dZ^ and, second, 
by inverting the Lagrange matrix to obtain the Pois- 
son matrix J = uj~^ with antisymmetric components 
jab ^ {Z", Z''}. Hence, we obtain the extended Poisson 
bracket defined in terms of two arbitrary functions F and 
G as {F, G} = daFr^dbG. 

The extended Vlasov equation is expressed in terms of 
the extended Vlasov distribution T{Z) and the extended 
Hamilton's equations as 







dT 
d7 



dZ" dT{Z) 



n}. 



(3) 



dr dZ^ 

In order to satisfy the physical constraint the ex- 
tended Vlasov distribution is required to be of the form 



T{Z) = cS[w-Hiz,t)] /(z,t). 



(4) 



where f{z,t) denotes the time-dependent Vlasov distri- 
bution on regular phase space. By integrating the ex- 
tended Vlasov equation (jSJ over the energy coordinate 
w (and using dr = dt), we obtain the regular Vlasov 
equation 







4f 
dt 



dl 
dt 



{/, H}. 



(5) 



Note that we use the same symbol { , } for the Pois- 
son bracket on regular phase space in Eq. © since w- 
derivatives appearing in the extended Poisson bracket 
vanish identically on regular phase space. 

Next, the extended Vlasov equation © is coupled with 
Maxwell's equations for the self-consistent electromag- 
netic fields 

IdE 47r 

V • E = 47r fl and V X B — = — J, (6) 
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where the charge-current densities 

- 1: (:) (7) 

are defined in terms of moments of the extended Vlasov 
distribution T (with d^p = c~^dw cfip) and the electric 
and magnetic fields E = — V$ — c~^dA/dt and B = 
V X A satisfy 

V-B = OandVxE + cr^dB/dt = 0. (8) 

The purpose of this paper is to show that: the asymptotic 
elimination of a fast time scale from the Vlasov equa- 
tion Q introduces polarization and magnetization effects 
into the Maxwell equations lO; and the reduced Vlasov- 
Maxwell equations possess exact conservation laws that 
can be derived from a variational principle by Noether 
method. 

The remainder of this paper is organized as follows. In 
Sec. ^ we present a brief summary of the Lie-transform 
perturbation method used in the asymptotic elimination 
of a fast time scale from the Vlasov equation and its asso- 
ciated Hamiltonian dynamics. In Sec. lIIII we present the 
reduced Vlasov-Maxwell equations introduced by the dy- 
namical reduction associated with a near-identity phase- 
space transformation. Through the use of the push- 
forward representation of particle fluid moments, we 
present expressions for the charge-current densities in- 
volving momentum-space moments of the reduced Vlasov 
distribution and present explicit expressions for the re- 
duced polarization charge-current densities and the di- 
vergenceless reduced magnetization current density. In 
Sec. IIVI we present the variational derivation of the re- 
duced Vlasov-Maxwell equations and derive, through the 
Noether method, the exact reduced energy-momentum 
conservation laws (explicit proofs are presented in Ap- 
pendix 0. Lastly, in Sec. we summarize the work 
present here and briefly discuss applications. 



II. DYNAMICAL REDUCTION BY 
NEAR-IDENTITY PHASE-SPACE 
TRANSFORMATION 

A. Near-identity Phase-space Transformation 

The process by which a fast time scale is removed from 
Hamilton's equations (^3) involves a near-identity trans- 
formation : Z —>■ Z{Z; e) = T^Z on extended particle 
phase space, where 

Z\Z,e) = Z^^ + eGl + e'I^G^, + (9) 



and its inverse near-identity transformation ^ : Z ^ 
Z(Z;e) = %'^Z, where 

Z'^(Z,e)^T-eGl~e'^G^,^^^^+---. m 

In Eqs. (j!!))- (|10|l . the dimensionless ordering parameter 
e ^ 1 is defined as the ratio of the fast time scale over a 
slow time scale of interest, and the nth-order generating 
vector field G„ is chosen to remove the fast time scale 
at order e" from the perturbed Hamiltonian dynamics. 
Examples of asymptotic elimination of fast time scales 
by Lie-transform perturbation method include guiding- 
center Hamiltonian theory gyrocenter Hamiltonian 
theory and oscillation-center Hamiltonian theory |^. 

B. Pull-back and Push-forward Operators 

Next, we define the pull-back operator on scalar fields 
induced by the near-identity transformation 

T,: T ^ T = T,J, (11) 

i.e., the pull-back operator Tg transforms a scalar field J- 
on the phase space with coordinates Z into a scalar field 
J- on the phase space with coordinates Z: 

T{z) = j,y{z) = y{%z) = y(z). 

Using the inverse transformation H1(J|I . we also define the 
push-forward operator: 

J-^ : T ^ y = J^^T, (12) 

i.e., the push-forward operator transforms a scalar 
field J- on the phase space with coordinates Z into a 
scalar field T on the phase space with coordinates Z: 

y(Z) ^ T-^T(Z) = T{T-^Z) = T{Z). 

Note that both induced transformations Hll() and H12|l 
satisfy the scalar- invariance property T{Z) = T(Z). 

In Lie-transform perturbation theory 3] , the pull-back 
and push-forward operators Hll|) and H12|l are expressed 
as Lie transforms: Tf^ = exp(±^^^-^ e" £„) defined in 
terms of the Lie derivative £„ generated by the nth-order 
vector field G„, which appear in the nth-order terms 
found in the near-identities lO and (|10|) . A Lie derivative 
is a special differential operator that preserves the tenso- 
rial nature of the object it operates on Q. For example, 
the Lie derivative £„ of the scalar field H is defined as 
the scalar field £„7i = G° daTi-, while the Lie derivative 
£n of a one-form F = Fa dZ° is defined as the one-form 

£n^ = G^ UJab^Z + d(G'", Fa), 

where Uab = dgTb ~ dbTa are the components of the two- 
form a; = dr. 

The pull-back and push- forward operators Hll|l - H12|) 
can now be used to transform an arbitrary operator 
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C : F{Z) — > C[F]{Z) acting on the extended Vlasov dis- 
tribution function T. First, since C[J^](Z) is a scalar 
field, it transforms to T7^{C[J^]}(Z) with the help of the 
push- forward operator 112|l . Next, we replace the ex- 
tended Vlasov distribution function T with its pull-back 
representation J- — T^J- and, thus, we define the trans- 
formed operator 



(13) 



By applying the induced transformation 113|l on the ex- 
tended Vlasov operator d/dr defined in Eq. , we obtain 



deF 
~d^ 



= t; 



dT 



T,F 



{F, nh 



(14) 



where the total derivative d^/dr along the transformed 
particle orbit is defined in terms of the transformed 
Hamiltonian 



and the transformed Poisson bracket 

{F, G}, EE T7i({T,F, T,G}). 



(15) 



(16) 



The Poisson-bracket transformation { , } — > { , }e can 
also be performed through the transformation of the ex- 
tended phase-space Lagrangian, = T^^F -I- dS, is ex- 
pressed as 



F, = Fo + e(Fi + dSi - £iFo) 



+ F2 + dS2 - £2^0 - £iFi 



1 



where S = e Si + e S2 + ■ ■ • denotes a (canonical) scalar 
field used to simplify the transformed phase-space La- 
grangian Fg at each order e" in the perturbation analy- 
sis. Note that the choice of S has no impact on the new 
Poisson-bracket structure 



dF, 



d (T-F) 



since d^5 — (i.e., d^f^S — d^^S — 0) and T^^ com- 
mutes with d. By inverting the reduced Lagrange ma- 
trix ZUe = 'T~^uj ^ Je = wj^, we, thus, obtain the 
reduced Poisson matrix J^, with antisymmetric compo- 
nents J° = Z }e, and define the reduced Poisson 

bracket {F, G}, = daF J, dbG. Lastly, we note that 
the extended-Hamiltonian transformation IjlSI) may be 
re-expressed in terms of the regular Hamiltonians H and 
H = J-^H - dS/dt as 



H 



Ha + £ { Hi — £iHq 



e H2 — £2Ho — £iHi 



_ dSi 
£\Hq 



dS2 
dt 



The new extended phase-space coordinates are chosen 
(i.e., the generating vector field G„ and the scalar field iS„ 
are specified at each order n = 1, 2, ... in the perturbation 
analysis) so that d^Z /dr — {Z , Ti.}e are independent 
of the fast time scale. 



III. REDUCED VLASOV-MAXWELL 
EQUATIONS 

A. Reduced Vlasov Equation 

The push-forward transformation of the extended 
Vlasov distribution yields the reduced extended 
Vlasov distribution 



T{Z) = cS[w-H{z,t)] /(z,t). 



(17) 



where the reduced extended Hamiltonian Ti, = H{z, t) — 
w is defined in Eq. (|15() . The extended reduced Vlasov 
equation 



d^T 



{T, n}e = 



(18) 



can be converted into the regular reduced Vlasov equa- 
tion by integrating it over the reduced energy coordinate 
w, which yields the reduced Vlasov equation 







dj 
dt 



dl 
dt 



{/, H}, 



(19) 



where /(z, t) denotes the time-dependent reduced Vlasov 
distribution on the new reduced phase space. Hence, we 
see that the pull-back and push-forward operators play 
a fundamental role in the transformation of the Vlasov 
equation to the reduced Vlasov equation. 



B. Reduced Maxwell Equations 

We now investigate how the pull-back and push- 
forward operators l|ll|l and 112|l are used in the transfor- 
mation of Maxwell's equations ©. The charge-current 
densities ((TJ can be expressed in terms of the general ex- 
pression (where time dependence is omitted for clarity) 



J^(r) = ^ e J d'^x J d^pv^'6'^{■K-v)T, 



(20) 



where the delta function (5'^(x — r) means that only parti- 
cles whose positions x coincide with the field position 
r contribute to the moment J'^(r). By applying the 
extended (time-dependent) phase-space transformation 
7^ : Z — > Z (where time t itself is unaffected) on the 
right side of Eq. H2UI) , we obtain the push-forward repre- 
sentation for J^: 

J^(r) = ^e j d^xd^p (T^^w'^) (5^ (x + - r) 



E 



d'^p 



(21) 



where ^v^ = (c, ^v) denotes the push-forward of 
the particle four-velocity v'^ and the displacement = 



4 



— X between the push- forward ^x of the particle 
position x and the (new) reduced position x is expressed 



(22) 



where the quadrupole contribution ^'S/ • [p^p^T) (which 
will be useful in what follows) is retained in Eq. 128|l . 
while the reduced electric-dipole moment (for each par- 
ticle species) 



e 



(29) 



in terms of the generating vector fields (Gi, G2, • • •) asso- 
ciated with the near-identity transformation. 

The push-forward representation for the charge- 
current densities, therefore, naturally introduces polar- 
ization and magnetization effects into the Maxwell equa- 
tions. Hence, the microscopic Maxwell's equations © 
are transformed into the macroscopic (reduced) Maxwell 
equations 



V-D = 47rp, 

^ 1 aD 47r - 

VxH — = —J, 
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(23) 
(24) 



where the reduced charge-current densities J = (cp, J) 
are defined as moments of the reduced Vlasov distribu- 
tion T: 



(cp, J) 



E 



^P-^ -dT 



(25) 



where d^^/dt denotes the reduced (e.g., guiding-center) 
velocity. The microscopic electric and magnetic fields E 
and B are, thus, replaced by the macroscopic fields 



D = E + 47rPe 
H = B 47rM, 



(26) 



where P,; and Me denote the polarization and magneti- 
zation vectors associated with the dynamical reduction 
introduced by the phase-space transformation 



C. Push-forward Representation of 
Charge-Current Densities 

We now derive explicit expressions for the reduced po- 
larization and the reduced magnetization by using 
the push-forward representation method. First, we de- 
rive the push-forward representation ()21f) for the charge 
density (J° = cp): 



p = ^ e j d'^pT 
= p - V-P„ 



E 



e I d p p^ T 



(27) 



where P = ^ / denotes the reduced charge den- 

sity that appears in Eq. H23(l and the polarization vector 
is defined as 



E 



d^p 



Pe - ^ V • ( p, p, 



(28) 



is associated with the fast-time-averaged charge separa- 
tion induced by the near-identity phase-space transfor- 
mation. 

Secondly, we derive the push-forward expression for the 
current density J, where the push- forward of the particle 
velocity v = dyi/dt (using the Lagrangian representation) 



T-i — 

' dt 



djX 

~(it 



^ dt ' 



t: 



dt 



(30) 



is expressed in terms of the reduced velocity d^li/dt, 
which is independent of the fast time scale, and the parti- 
cle polarization velocity d^p^/dt, which has both fast and 
slow time dependence. Note that the fast-time-average 
particle polarization velocity d^p^/dt, which is nonvan- 
ishing under certain conditions, represents additional re- 
duced dynamical effects (e.g., the standard polarization 
drift in guiding-center theory 0, ^2|) not included in 
djc/dt. Hence, the push-forward expression 1)21(1 for the 
current density J is 



dt J 



•C31) 



We may now replace the polarization velocity d^p^/dt 
in Eq. H31|l by using the following identity based on the 
reduced polarization vector if^ : 



dt ^ 



E^ 



4pA y 



dt 
d'^pT 



dfX 
~dt 



(32) 



where the reduced Vlasov equation (|14|l was used and in- 
tegration by parts was performed. Using the vector iden- 
tity V • (BA - AB) = V X (A X B), the push-forward 
representation for the current density is, therefore, ex- 
pressed as 



9Pe 

dt 



c V X 



(33) 



where J = ^ / d'^P (d^x/dt) T denotes the reduced 
current density appearing in Eq. (Pljl . Jpoi = dV^^/dt 
denotes the reduced polarization current, and the di- 
vergenceless reduced magnetization current Jmag = 



5 



c V X Me is expressed in terms of the reduced magne- 
tization vector 

+ (34, 

which represents the sum (for each particle species) of the 
intrinsic (fast-time-averaged) magnetic-dipole moment 



e / 



2c 



dt 



(35) 



and a moving electric-dipole contribution (tt^ X dfxjdt), 
as suggested by classical electromagnetic theory 



0. 



IV. VARIATIONAL FORMULATION OF 
REDUCED VLASOV-MAXWELL EQUATIONS 



and we introduce the reduced antisymmetric 
polarization-magnetization tensor |l7j 



4.^ 9^ 



(39) 



where the reduced polarization and magnetization vec- 
tors = PI and K'^ = e'^'' k are defined as 

(P.,M,,.-i:/..p7(§||). ,40, 

We begin with an expression for the variation of the 
reduced Lagrangian density 



47r 



^ J d^pS [t, h} 



(41) 



d 



An — 1/ 
c 



A. Reduced Variational Principle 

We now show that the reduced Vlasov-Maxwell equa- 
tions iCHll and (ESl-lEl 

can be derived from a varia- 
tional principle J d'^x 5C = 0. Variational principles for 
reduced Vlasov-Maxwell equations have been presented 
previously by Pfirsch Ts] and Pfirsch and Morrison 
using the Hamilton-Jacobi formulation and by Kaufman 
et al. 0, 0, 0, 0, 01 using the Low-Lagrangian for- 
malism. Here, we present a variational principle for re- 
duced Vlasov-Maxwell equations based on the reduced 
Lagrangian density |23| 



F : F 



J2 / d^P^K 



(36) 



where and TC denote the reduced extended Vlasov 
distribution and the reduced extended Hamiltonian, re- 
spectively. In this Section, we use the convenient space- 
time metric g^'^ — diag (— 1, 1, 1, 1), so that the electro- 
magnetic field tensor F^i, = df^A^ — d^Af^ is defined in 
terms of the four-potential A^^ = (— $, A) and, thus, 
Fio = Ei and Fij = e^fc . In order to simplify our pre- 
sentation, we use canonical four-momentum coordinates 
= (— UJ/c, p), so that the reduced extended Hamilto- 
nian Ti. is required to be invariant under the gauge trans- 
formation 



A„ 



+ d^,x and p 



+ {elc)d^x, (37) 



where x denotes the electromagnetic gauge field. 

Note that, as a result of the dynamical reduction of 
the Vlasov equation, the reduced Hamiltonian TL is not 
only a function of the four-potential A^ but also of the 
field tensor F^^. From these dependences, we express the 
reduced four-current density (I25|l as 



r 



(cp, J) 



dn 



(38) 



which is generated by the four-potential variation SA^, 
and the Eulerian variation for the reduced Vlasov 
distribution ST = — 5Z daT = {S, J-}^, where S is 
the generating scalar field for a virtual displacement on 
reduced phase space, 5Z = {Z , S}^. Note that the 
divergence term d^J^ , where the reduced Noether four- 
density is 



deX^ SAy 



dt 



in 



AnK''^' 



does not contribute to the variational principle 
/ d'^x SC = but instead is used to derive exact con- 
servation laws by applications of the Noether method. 

Next, as a result of the variational principle, where 
the variations S and SA^, are arbitrary (but are required 
to vanish on the integration boundaries), we obtain the 
reduced Vlasov equation lfTH|l . {T, Ti}e = 0, and the 
reduced (macroscopic) Maxwell equations 



dxt^ ^ ' c 



(43) 



from which we recover the reduced Maxwell equations 
(|23|l and (|24() . Here, the polarization-magnetization four- 
current is expressed in terms of the tensor (|39|l as 

-^^,K^''' - (-V-Pe,c-iatPe + VxM,) 

— (Ppob C Jpol + C Jmag) ■ 

Note that the electromagnetic field tensor also satisfies 
daF^i, + df_iFi,a + di,F„fj_ — 0. In addition, we note that 
the reduced electric-dipole and reduced magnetic-dipole 
moments (|29|) and H35|) are also expressed in terms of 
derivatives of the reduced Hamiltonian H as 



dH 



dH 1 4x dH'' 
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which provides a useful consistency check on the reduced 
Vlasov-Maxwell equations. Lastly, the reduced charge 
conservation law 

follows immediately from the reduced Maxwell equations 
(|43|1 as a result of the antisymmetry of i^^" and K^'^ (i.e., 

B. Reduced Energy- momentum Conservation Laws 



Here, the antisymmetric tensor K^" , defined in Eq. H39|l . 
represents the effects of reduced polarization and mag- 
netization. Note that, while the last two terms are in- 
dividually gauge-dependent, their sum is invariant under 
the gauge transformation H37I) . Explicit proofs of energy- 
momentum conservation for the reduced Vlasov-Maxwell 
equations based on the reduced energy-momentum ten- 
sor (|4ti|) are presented in Appendix ^ Lastly, addi- 
tional angular-momentum conservation laws can be de- 
rived from the reduced Noether equation (5£ = d^^J^ by 
considering invariance of the reduced Lagrangian density 
L, with respect to arbitrary rotations in space. 



We now derive the energy-momentum conservation law 
from the reduced Noether equation (5£ = d^i^J^ associ- 
ated with space-time translations generated by bx" = 
(ci5<, (5x), where the variations {S,6A^^SC) are 



SUMMARY 



SA„ 



SjC 



d,{A„Sx'') ) 



da (jC Sx" 



(45) 



Here, we note that 5x^ = {x'^, S}e, while the variations 
6A1, and 6C are expressed in terms of the Lie derivative 
£sx generated by Sx" as SA^^dx'^ = — £sx{A,y dx^) and 
S£d^Z = - £sx(Cd^Z). 

After some cancellations introduced through the use 
of the reduced Maxwell equations H43() . we obtain the 
reduced energy-momentum conservation law d^T^^'^ = 0, 
where the reduced energy-momentum tensor is defined as 



16tt 



F : F 



An ^ 



47r K^"')F„ 



dt 



In this paper, the general theory for the reduced 
Vlasov-Maxwell equations was presented based on the 
asymptotic elimination of fast time scales by Lie- 
transform Hamiltonian perturbation method. This dy- 
namical reduction is based on a near-identity transfor- 
mation on extended phase space, which induces transfor- 
mations on the Vlasov distribution and the Vlasov op- 
erator, as well as introducing a natural (push-forward) 
representation of charge-current densities in terms of re- 
duced charge-current densities and their associated re- 
duced polarization and magnetization effects. The varia- 
tional formulation of the reduced Vlasov-Maxwell equa- 
tions allows the derivation of exact energy-momentum 
conservation laws by Noether method. 

The Table shown below summarizes the polarization 
and magnetization effects observed in reduced Vlasov- 
Maxwell equations that have important applications in 
plasma physics. 



Reduced Dynamics 






Guiding-center 




— Jib 


Gyrocenter 


e (cbo/Bo) X (eAi/c + musi + P|| Bi/Bq) 


-7Z(bo + eBi/So) 


Oscillation-center 


e^ek X (-i|* X |) 


e^euj'/c {-if x|) 



r 



First, in guiding-center Hamiltonian theory 0, |^ background electric field E = — V$, the fast time scale 
for a strongly magnetized plasma in the presence of a is associated with the rapid gyromotion of a charged par- 
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tide about a magnetic field line and fast-time-averaging 
is carried out by averaging with respect to the gyroan- 
gle. The fast-time-averaged reduced electric-dipole mo- 
ment Tfgc = epg(, includes effects due to the background 
electric field (ve — cE X 'B/B'^) as well as the magnetic 
(VB and curvature) drift velocity vb- The fast-time- 
averaged particle polarization velocity dgcVgc/ dt includes 
the standard polarization drift velocity {c/ Bi^) d'E±/dt, 
which is not included in the guiding-center drift veloc- 
ity dgcx/dt = Vgc — Ve -l- Vb. On the other hand, the 
reduced magnetic-dipolc moment yields the classical par- 
allel magnetization term /Jgc = — Jib (where ]I denotes 
the guiding-center magnetic-moment adiabatic invariant 
and b = Ji/B denotes the unit vector along a magnetic- 
field line), which enables the reconciliation of the parti- 
cle current J with the guiding-center current Jgc through 
the relation J = Jgc -I- Jmag (valid in a static magnetized 
plasma) . 

Next, gyrocenter Hamiltonian theory describes 
the reduced (gyroangle-independent) perturbed guiding- 
center Hamiltonian dynamics associated with low- 
frequency, electric and magnetic fluctuations (eEi,eBi) 
in a strongly-magnetized plasma (with static magnetic 
field Bo = -Bo bo and Eo = 0). Note that the re- 
sults shown here are valid only in the limit of zero 
Larmor radius 0- The reduced electric-dipole mo- 
ment includes not only the perturbed polarization-drift 
term (mc bo/Bo X uei), but also the effects due to 
magnetic flutter (cp|| bo X Bi/i?o) and the inductive 
part of the perturbed E x B velocity (i.e., the polar- 
ization drift velocity includes the higher-order correc- 
tion —dtAi X bo/i?o to the perturbed E x B veloc- 
ity — V$i X cbo/i?o). On the other hand, the reduced 
magnetic-dipole moment includes a correction to the clas- 
sical parallel magnetization term due to the perturbed 
magnetic field: Jl^y = — JI {bo + eBi/ Bq), where Jl now 
denotes the gyrocenter magnetic-moment adiabatic in- 
variant. 

Lastly, oscillation-center Hamiltonian theory Q de- 
scribes the reduced dynamics of charged particles inter- 
acting with a high-frequency electromagnetic wave in a 
weakly-inhomogeneous plasma for which the eikonal ap- 
proximation is valid. The eikonal representation for the 
wave fields is (Ei,Bi) = (Ei,Bi) exp(zeQ^O) -I- c.c, 
where eo 1 denotes the eikonal small parameter 
while the eikonal phase 8(eor, coO is used to define 
uj = —e^^dfO and k = Cg^^VO, with uj' = w — k-v 
denoting the Doppler-shifted wave frequency. Note that 
fast-time-averaging, here, is carried out by averaging with 
respect to the eikonal phase 0. By considering the sim- 
plest case of an unmagnetized plasma 0, the first-order 
term for the displacement Pe = e ^ + • • • has the eikonal 
amplitude 

V 



El + - xBi , 



where e <^ 1 denotes the amplitude of the electromag- 



netic wave represented by the first-order fields Ei and Bi. 
Note that both reduced oscillation-center electric-dipole 
and magnetic-dipole moments are quadratic functions of 
the wave fields. An additional wave-action conservation 
law results from the invariance of the reduced Lagrangian 
density on the eikonal phase but its derivation is outside 
the scope of this work. 
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APPENDIX A: EXPLICIT PROOFS OF 
ENERGY-MOMENTUM CONSERVATION FOR 
THE REDUCED VLASOV-MAXWELL 
EQUATIONS 

In this Appendix, we present explicit proofs of energy- 
momentum conservation based on the reduced energy- 
momentum stress tensor (|46|) . We begin with the reduced 
energy conservation law d£/dt -I- V • S = 0, where the 
reduced energy density £ = T^" is 

f = - ^ (|Ep - IBp) + ^ 

where the w-integration was performed. First, using the 
reduced charge conservation law (|44|l and the identity 



- dH 

dt 



a$ J dA 



dt 



dt 



^ dE 9B 



we obtain 

d£ _ E dD U dB 
'dt ~ i^'~dt ^ i^'~dt 



E- J 



V- (J$) - d'pH{f, H}, 



where we used the reduced Vlasov equation ((T^ to ob- 
tain the last term. Lastly, using the reduced Maxwell 
equations (|SJ) and 124|) and the identity 



d-'pHif, H} - / d^p{fH, H} 



= V. 



d'pfH 



^d. 



dt 



we finally obtain 

d£_ 
'dt 



V- f— ExH - J$ 



8 



where the energy-density flux S*' = cT*". 

Next, we consider the reduced momentum conservation 
law dH/dt + V • T = 0, where the reduced momentum 
density W = T'^^/c is 



n 



D X B 

An c 



which has the Minkowski form T^. By substituting the 
reduced Vlasov-Maxwell equations © and H23|) - (|24|l and 
the reduced charge conservation law H44() . we obtain 



an 

'dt 



(BH 



DE) 



4tt 

(VH-B + VE-D 

Lastly, using the identities 



JA 



V$p - VA-- 

c 



and 



VHB + VED 



we finaly obtain 



dn 

'dt 



VE-P, - VB-Me, 



lEl' 



iBr 



47rM, -B 



2 2 
47r(VE-P, + VB-M,), 



1 



An 

+ - JA^ 

c An 

- V-T. 



lEl 



IBI 



^ dt 
47r M, • B 



^ j d^p fVH = V$ p - VA ■ 
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